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Abstract 

The maize (Zea mays L.) Aux/IAA protein RUM1 (ROOTLESS WITH UNDETECTABLE MERISTEMS 1) controls seminal 
and lateral root initiation. To identify RUM1 -dependent gene expression patterns, RNA-Seq of the differentiation zone 
of primary roots of ruml mutants and the wild type was performed in four biological replicates. In total, 2 801 high- 
confidence maize genes displayed differential gene expression with Fc >2 and FDR <1%. The auxin signalling-related 
genes rum1, Iike-auxin1 [Iax1), Iax2, [nam atafcuc 1 nad), the plethora genes plt1 {plethora 1), bbml (baby boom 1), 
and hscfl {heat shock complementing factor 1) and the auxin response factors arf8 and arf37 were down-regulated in 
the mutant ruml. All of these genes except nad were auxin-inducible. The maize arf8 and arf37 genes are orthologues 
of Arabidopsis MP/ARF5 (MONOPTEROSIARFS), which controls the differentiation of vascular cells. Histological anal- 
yses of mutant ruml roots revealed defects in xylem organization and the differentiation of pith cells around the 
xylem. Moreover, histochemical staining of enlarged pith cells surrounding late metaxylem elements demonstrated 
that their thickened cell walls displayed excessive lignin deposition. In line with this phenotype, rum 1 -dependent 
mis-expression of several lignin biosynthesis genes was observed. In summary, RNA-Seq of RUM1 -dependent gene 
expression in maize primary roots, in combination with histological and histochemical analyses, revealed the specific 
regulation of auxin signal transduction components by RUM1 and novel functions of RUM1 in vascular development. 

Keywords: Auxin, lateral roots, lignification, RNA-Seq, RUM1 , vasculature, xylem. 



Introduction 

Maize (Zea mays L.) plays an important agronomic role 
as feed, food, and source of bioethanol. The complex root 
system of maize facilitates water and nutrient uptake and 
anchorage of the plant (Aiken and Smucker, 1996). The 
embryonic root system of maize consists of a primary root 
and a variable number of seminal roots. The post-embry- 
onic root system comprises lateral and shoot-borne roots. In 
all maize root types lateral roots are initiated from phloem 
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pole pericycle and endodermis cells (Hochholdinger et al, 
2004a). 

The plant hormone auxin is a key regulator of many aspects 
of plant development. In roots, auxin controls lateral root 
development (Jansen et al, 2012; Peret et ah, 2009a, b), root 
architecture (Overvoorde et al, 2010) and vascular develop- 
ment (Wilson et al, 1991; Uggla et al, 1996; Mattsson et al, 
1999). It was suggested that polar auxin transport controls 
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lateral root initiation in maize (Jansen et al, 2012). Moreover, 
an auxin maximum was observed from late metaxylem to pro- 
toxylem in cells surrounding protophloem cells in the trans- 
genic maize marker line DR5-RFP (Jansen et al, 2012). This 
auxin maximum is likely to be essential for lateral root position- 
ing and initiation. By contrast, inhibition of polar transport 
prevents an auxin response maximum from xylem to phloem 
pole cells in an anticlinal orientation and results in disorgan- 
ized vascular tissues (Mattsson et al, 1999; Jansen et al. , 2012). 

Central regulators of auxin signalling include the Transport 
inhibitor -Response 1 (TIR1) protein, members of the TIR1- 
like ^4uxin F-Box (AFBs) family, ^4M.vin//ndole Acetic ^4cid 
(Aux/IAA) proteins, and ^4uxin Response factor (ARF) pro- 
teins (Mockaitis and Estelle, 2008). At low intracellular auxin 
concentrations, Aux/IAA proteins act as transcriptional 
repressors that interact with ARF proteins via their domains 
III and IV. The ARF proteins of these complexes interact with 
Auxin Responsive Elements (AuxREs) in the promoters of 
downstream genes, thereby repressing their transcription. By 
contrast, high auxin levels stabilize interactions between Aux/ 
IAA proteins and SCF TIR1 E3 ubiquitin-ligase complexes, 
resulting in the degradation of Aux/IAA proteins by the 26S 
proteasome (Gray et al, 2001; Tian et al, 2002; Woodward 
and Bartel, 2005; Tan et al, 2007). As a consequence, ARF 
proteins released from the Aux/IAA interactions can promote 
auxin-responsive target gene transcription. 

The semi-dominant maize ruml {rootless with undetectable 
meristem 1) mutant is blocked in the initiation of embryonic 
seminal roots and post-embryonic lateral roots of the primary 
root (Woll et al, 2005). The mutant ruml displays several 
auxin-related defects in root development. While exogenous 
auxin (aNAA) application induced additional lateral roots 
in wild-type primary roots, aNAA did not initiate any lateral 
roots in ruml (Woll et al, 2005). Moreover, polar auxin trans- 
port in ruml primary roots was reduced by 83% compared 
with wild-type primary roots while it was not affected in the 
coleoptile of the mutant ruml (Woll et al, 2005). The ruml 
gene (Genbank AC: GRMZM2G037368) encodes ZmlAAlO, 
a member of the Aux/IAA protein family (Wang et al. ,2010). 
RUM1 interacts with the transcriptional activators ARF25 
and ARF34 (von Behrens et al. , 201 1 ). The mutated ruml pro- 
tein lacks 24 amino acids including the degron motif 'GWPPV 
in domain II of ruml. Therefore, ruml cannot interact with 
the SCF TIR1 E3 ubiquitin-ligase complexes which prevents its 
ubiquitin-mediated proteasomal degradation and resulting in 
constitutive repression of downstream gene expression (von 
Behrens et al, 2011). In the present study, genes differentially 
expressed between wild-type and ruml primary roots were 
identified via RNA-Seq suggesting direct or indirect regula- 
tion of these genes by RUML In combination with histo- 
logical and histochemical analyses, a RUM 1 -dependent gene 
network was identified and novel functions of RUM1 in vas- 
cular development were revealed. 

Materials and methods 

Plant material and growth conditions 

Seeds of the F 8 -generation of the maize mutant ruml and its homozy- 
gous wild type obtained by seven cycles of selling of heterozygous 



plants were used in these experiments. Seeds were sterilized with 6% 
hypochlorite under vacuum at 500 mPa for 5 min, rinsed five times 
in distilled water, and germinated in paper rolls in a plant growth 
chamber at 28 °C with a 16/8 h light/dark regime at 21 °C (Woll et al. , 
2005). For the RNA-Seq experiment, primary roots, 2cm in length, 
of the mutant ruml and its homozygous wild type were harvested. 
Subsequently, 5 mm of the root tip, including the meristematic and 
elongation zones, were removed. Hence, only the differentiation 
zone of these roots was analysed by RNA-Seq. For auxin induc- 
tion experiments, 5-d-old seedlings of the maize inbred line B73 that 
were grown under the same conditions were treated with the auxin 
analogue aNAA (a-naphthyl acetic acid; working solution 5 uM) 
for 3 h. The differentiation zone of the ~5 cm primary roots was har- 
vested after 0, 1, 2, and 3h of aNAA exposure (Taramino et al, 
2007) for subsequent analyses. For histological experiments, wild- 
type and ruml seedlings were grown in aNAA (0. 1 uM) or the auxin 
transport inhibitor NPA (1-iV-naphthylphthalamic acid, 10 uM). 

RNA extraction and RNA sequencing 

Total RNA was extracted via the RNeasy Plant Mini Kit (Qiagen, 
Hilden, Germany) and 3^4 roots were pooled per replicate, sub- 
sequently treated with RNase-free DNasel (Fermentas, St Leon- 
Roth, Germany). RNA quality of all samples was assayed via a 
Bioanalyzer (Agilent Technologies, Boeblingen, Germany). As sug- 
gested by Agilent only samples with an RIN (RNA integrity number 
>7 were used for downstream analyses. Per genotype (wild type ver- 
sus ruml) or treatment (auxin induction) four biological replicates 
were surveyed. The cDNA libraries for RNA-Seq were constructed 
using the TruSeq RNA sample preparation kit (Illumina Inc., San 
Diego, CA, USA). Each of the four replicates of each genotype was 
bar coded with one of the Illumina indices AR001, AR003, AR008, 
and AR009. Subsequently, the four wild-type replicates were pooled 
and loaded onto lane three of a flow cell, while the four ruml repli- 
cates were pooled on lane four of a flow cell. These lanes were then 
sequenced using a Genome Analyzer II (Illumina Inc., San Diego, 
CA, USA) according to the manufacturer's instructions resulting in 
146 bp single-end reads. 

RNA-Seq mapping and statistical analysis 

RNA-Seq reads were trimmed according to Liu et al. (2012) and sub- 
sequently mapped to the B73 reference genome (ZmB73_RefGen_v2) 
(Schnable et al., 2009) using GSNAP (Wu and Nacu, 2010). For sub- 
sequent analyses, only reads mapping uniquely (<2 mismatches every 
36 bp and fewer than five bases for every 75 bp as tails) to the refer- 
ence genome were used. Genes with an average of at least five uniquely 
mapped reads across libraries and at least two samples with positive read 
counts («=22 833) were tested for differential expression using the R 
package QuasiSeq (http://cran.r-project.org/web/packages/QuasiSeq). 
The 0.75 quartile of reads from each library was used as the normali- 
zation factor (Bullard et al, 2010). The negative binomial QLSpline 
method implemented in the QuasiSeq package was used to compute 
a p- value for each gene (Lund et al, 2012). A multiple test controlling 
approach (Nettleton et al. , 2006) was used to convert the p- values to 
q- values (Storey, 2002). Genes with FDR <1% and Fc >2 were declared 
to be differentially expressed (Li et al, 201 1). Subsequently, genes were 
classified into functional categories via MapMan (Thimm et al. , 2004). 
To determine if specific functional groups are overrepresented among 
the differentially expressed genes with reference to all expressed genes, 
the expected number of genes for each functional category was cal- 
culated based on the distribution of functional categories among all 
expressed genes. To determine if significantly more or less genes than 
expected were detected for each individual category a % 2 test for inde- 
pendence with Yates' continuity correction was performed. 

qRT-PCR expression analyses 

cDNA for qRT-PCR was synthesized from 500 ng total RNA 
using the qScript cDNA SuperMix (Quanta Biosciences, 
Gaithersburg, MD, USA). qPCR was performed in a CFX384 
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Real-Time PCR detection system (Bio-Rad, Munich, Germany) 
for each of the four biological replicates in three technical rep- 
lications in a total reaction volume of 8 pi using the MESA 
Green qPCR Mastermix Plus for SYBR Assay no ROX kit 
(Eurogentech, Cologne, Germany). Primers with an efficiency 
between 0.9 and 1.1 were used for qPCR, which was tested in 
a dilution series (1, 1/2, 1/4, 1/8, 1/16, 1/32, 1/64, 1/128). Gene 
expression for each genotype and each time point of auxin induc- 
tion was assayed relative to myosin (Genbank AC: 486090G09. 
xl; primers: 486090G09.xl-5'; 486090G09.xl-3') which has been 
previously used as a qPCR standard for expression analyses in 
maize roots (von Behrens et at, 2011). The oligonucleotide prim- 
ers ruml-fw and ruml-rv (ruml, GRMZM2G037368), ar/S-fw 
and arf8-rv (arJS, GRMZM2G034840), arf37-fw and ar/37-rv 
(arf3 7, GRMZM2G086949), nacl -fwand«ac7 -r\(nacl , GRMZM2 
G081930),/ax7-fwand/ajc7-rv(/fl.\-7,GRMZM2G129413),/ax2-fw 
and lax2-rv (lax2, GRMZM2G149481), pltl-fv/ and pltl-rv (phi, 
GRMZM2G141638), bbml-fw and bbml-rv (bbml, GRMZM2 
G366434), hscfl-iw and hscfl-rv (hscfl, GRMZM2G139082), 
cad-fw and cad-xv (cad, GRMZM2G443445), and /5/z-fw and 
f5h-r\ (f5h, GRMZM2G100158) were used for surveying expres- 
sion of these genes (see Supplementary Table S5 at JXB online). 
Differential gene expression was determined by Student's / tests 
(*p <0.05; **p <0.01; ***p <0.001; «=4). 

Histology and histochemistry 

Feulgen staining of whole roots was performed as previously 
described by Woll et al. (2005). For histological analyses and subse- 
quent lignin staining, root fragments of the differentiation zone of 
the primary root were fixed in 4% paraformaldehyde in 1 x phosphate 
buffer (containing 0.01 M Na 2 HP0 4 and 0.01 M NaH 2 P0 4 , pH 7) for 
2h under vacuum at 100 mPa. Subsequently, the root fragments were 
sectioned by hand, mounted with water and then transferred into 
65% glycerol. Transverse sections were analysed under bright field 
using a dissection microscope (PixCell He System, Zeiss). For lignin 
staining, transverse sections were incubated in 1% phloroglucinol in 
12% HC1 for lOmin, then transferred into acid solution (75% glycerol 
in 10% H 2 S0 4 ) and observed under bright field conditions using a 
dissection microscope (PixCell He System, Zeiss). 

Results 

Comparative RNA-Seq analysis of ruml and wild-type 
primary roots 

The maize ruml mutant is impaired in lateral root initiation 
in the primary root. To understand the molecular network 
regulated by the Aux/IAA protein RUM1, ~2cm long primary 
roots of the mutant ruml and homozygous wild-type siblings 
were subjected to an RNA-Seq analysis (Fig. 1A). At this early 
developmental stage no lateral root primordia were detectable 
in primary roots of wild-type seedlings by Feulgen staining 
(data not shown). Hence, at this stage, primary roots of wild- 
type and ruml were morphologically indistinguishable. Lateral 
roots in maize are formed in the differentiation zone of roots. 
Therefore, only the differentiation zone of primary roots, i.e. 
primary roots, after the removal of 5 mm of the root tip which 
comprises the meristematic zone and most of the elongation 
zone, were used for subsequent transcriptome profiling. A flow 
chart of the RNA-Seq experiment is provided in Fig. IB. 
Sequencing of four biological replicates of wild-type and ruml 
libraries yielded, on average, ~1 5 million 146 bp single-end reads 
per sample (see Supplementary Table SI at JXB online). After 
removal of bar-code tags from ruml and wild-type reads, on 
average, 76% and 73%) of the trimmed reads mapped to unique 



positions of the maize B73 reference genome (ZmB73_RefGen_ 
v2; see Supplementary Table SI at JXB online), respectively. Of 
these reads, 94-95%o were associated with high-confidence gene 
models of the 'filtered gene set' (4a.53_v2; maizesequence.org). 

Identification of differentially expressed genes 

The RNA-Seq reads from the differentiation zone of 2-cm- 
long maize primary roots aligned to 22 833 of 39 656 
(58%o) genes in the high confidence 'filtered gene set' (see 
Supplementary Table S2 at JXB online). The distribution of 
fold-changes versus the mean of gene expression in the dif- 
ferentiation zone of primary roots between ruml and the 
wild type demonstrated that differential gene expression was 
observed for genes with both low and high expression levels 
(Fig. 1C). In total, 8 309 (FDR <5%>), 5 701 (FDR <1%), and 
3 502 (FDR <0.1%o) genes displayed differential expression at 
different significance levels irrespective of their fold-changes 
(Fig. ID). Hence, these classes represented 36%o, 25%o, and 
15%> of all genes for which expression was detected. At FDR 
<l%o, 2 801 (12% of all) genes displayed a Fc (fold change) >2 
(see Supplementary Table S2 at JXB online), including 1 741 
genes that were down-regulated (Fig. IE, red dots) and 1 060 
genes that were up-regulated in the ruml mutant (Fig. IE, blue 
dots). Of those genes, 1 723 (62%o) were functionally annotated 
using the MapMan software (see Supplementary Table S3 at 
JXB online), including 205 genes that were assigned to more 
than one functional class. To identify over- and underrepre- 
sented functional classes among the differentially expressed 
genes, the expected number of genes for each functional class 
was calculated based on all expressed genes and subsequently 
compared with the detected number of genes in this class. This 
analysis revealed that significantly more genes than expected 
(p <0.001) were assigned to the classes photosynthesis, cell 
wall, metal handling, secondary metabolism, miscellaneous, and 
DNA (see Supplementary Table S4 at JXB online). 

qRT-PCR confirmation of differentially expressed genes 

To confirm the RNA-Seq results independently, expression 
levels of a subset of differentially expressed genes related 
to auxin signalling were analysed by qRT-PCR. Consistent 
with prior results (von Behrens et al., 2011), the Aux/IAA 
gene ruml (GRMZM2G037368), used as a positive con- 
trol, was down-regulated in the ruml mutant compared with 
the wild type (Fig. 2A). Moreover, three plethora (pit) genes 
pltl (GRMZM2G141638), bbml (GRMZM2G366434), 
and hscfl (GRMZM2G139082) were down-regulated in 
ruml (Fig. 2A). Furthermore, two auxin response factors 
arf8 (GRMZM2G034840) and arf37 (GRMZM2G086949), 
two auxin influx transporters like-auxl (laxl, 
GRMZM2G129413) and lax2 (GRMZM2G149481), and 
a nacl (GRMZM2G081930) gene were down-regulated in 
ruml compared with wild-type primary roots. 

Auxin inducibility of differentially expressed genes 
involved in auxin signal transduction 

To test the auxin inducibility of these auxin signal trans- 
duction genes, time-course experiments were performed in 
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Fig. 1. (A) Phenotype of 2cm primary roots of homozygous wild type (left) and ruml (right). The differentiation zone of the primary roots was 
subjected to RNA-Seq experiments. Lateral root primordia do not develop at this stage (Woll ef al. , 2005). Scale bar: 1 cm. (B) Flow chart of the 
lllumina RNA-Seq experimental design. (C) MA-similar plot providing an overview of the correlation of gene expression changes versus expression 
intensity by plotting the log 2 of the fold change of rumJ/wild type against the log 2 of the mean of gene expression in ruml and the wild type from 
RNA-Seq data. Each dot represents an expressed gene. (D) Histogram of FDR-values resulting from the comparison of gene expression between the 
ruml mutant and wild-type primary roots. (E) Volcano plot displaying differential gene expression patterns between ruml and the wild type. The /-axis 
denotes negative log 10 FDR of each expressed gene. The X-axis denotes log 2 -fold changes of rum7/wild type. Each dot represents an expressed 
gene. Red dots denote down-regulated genes in ruml , blue dots denote up-regulated genes in ruml with cut-off thresholds of Fc >2; FDR <1%. (This 
figure is available in colour at JXB online.) 



Control of root vascular differentiation by RUM1 | 4923 



arf37 




plt1 




RNA-Seq: Fc = 49 




RNA-Seq: only WT 




RNA-Seq: only WT 



hscfl 



hC - 3Uf 

RNA-Seq: Fc = 76 



B 





2h 3h 




RNA-Seq: Fc = 13 



Fc 
* 

II 


= 1.3 
* 

h 

rih 


Oh 1h 2h 3h 

bbml 



r 1- 



2h 3h 



!ax1 




Oh 1h 2h 3h 



Oh 1h 2h 3h 




, 0.04 
0.03 
0.02 



arf37 



n 




1h 2h 



UJ. 




1h 2h 3h 



Fig. 2. (A) Confirmation of differential gene expression of nine genes involved in auxin signal transduction using qRT-PCR. Each of these genes was 
differentially expressed between rum1 and the wild type in the RNA-Seq experiment. In each graph fold-changes (Fc; rum7/wild type) from the qRT-PCR 
and the RNA-Seq experiment are displayed. (B) Auxin inducibility was tested for the nine auxin-related genes by a qRT-PCR time-course experiment with 
reference to time point 0 h. Wild-type primary roots (5-d-old) were treated with aNAA and harvested at 0, 1,2, and 3 h of exposure. Statistical analyses 
were performed with a two-sided Student's f test (*p <0.05; **p <0.01 ; ***p <0.001 ). 
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5-d-old primary roots treated with 5 uM aNAA (Fig. 2B). 
Expression was measured relative to the housekeeping gene 
myosin (Genbank AC: 486090G09.xl) for each time point. 
The genes ruml, laxl, and lax2 and all three pit genes (pltl, 
bbml, and hscfl) were induced within 3h of aNAA exposure. 
Remarkably, expression of arf8 and cirf37 was induced within 
2 h of aNAA exposure but returned to control levels 3 h after 
auxin treatment. The transcription of nacl was not affected 
by aNAA exposure. 

RUM1 regulates vascular development 

The homoeologues arf8 (GRMZM2G034840) and arJ37 
(GRMZM2G086949), which were significantly down-regu- 
lated in ruml (Fig. 2 A) encode putative transcriptional acti- 
vators. Their closest homologue in Arabidopsis is MPIARF5, 
which regulates vascular development (Przemeck et al, 1996; 



Hardtke and Berleth, 1998). To determine whether ruml 
mutants are impaired in vascular development, a comparative 
histological analysis of wild-type and ruml primary roots was 
performed at different developmental stages. In both geno- 
types, roots of 20, 40, 60, and 80mm length were dissected by 
serial transverse hand sections in 10 mm increments (Fig. 3). 
In wild-type primary roots, at all four developmental stages 
(20, 40, 60, and 80 mm), protoxylem elements, early metax- 
ylem elements, and parenchymous pith cells surrounding the 
late metaxylem elements were well developed and properly 
arranged throughout development (Fig. 3A-T). By contrast, 
in primary roots of the mutant ruml, xylem elements and 
the pith cells were largely in disarray. Only the oldest part of 
the mutant roots at all four developmental stages displayed 
properly arranged xylem elements (20 mm: Fig. 3A'; 0 mm: 
Fig. 3C; 60 mm: Fig. 3G'; 80 mm: Fig. 3M'). Moreover, in 
60 mm and 80 mm ruml roots, only the oldest parts of the 
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Fig. 3. Series of transverse sections displaying the central cylinder of primary roots at different developmental stages: 20 mm wild-type (A, B) and 
ruml (A', B') primary roots, 40mm wild-type (C-F) and ruml (C'-F') primary roots, 60mm wild-type (G-L) and ruml (G'-L') primary roots, 80mm wild- 
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primary root (60 mm: Fig. 3G-T; 80 mm: Fig. 3M -O') dis- 
played well differentiated xylem elements, while the younger 
parts of these roots (60 mm: Fig. 3K'-L'; 80 mm: Fig 3P / -T / ) 
did not show any differentiated xylem elements. Hence, the 
magnitude of alterations in the vasculature of ruml corre- 
lated with the developmental stage at which these parts of 
the root were formed. In all analysed mutant roots, vascu- 
lar differentiation was not affected in cells that were released 
from the meristematic zone soon after germination (Fig. 3 A', 
C, G', M') whereas the defects became more severe the later 
these cells were released from the meristem. 

The number of xylem elements and the number and size of 
pith cells was determined in representative cross-sections of 
80 mm ruml roots at a distance of 70 mm from the root tip. In 
these sections, the number of early and late metaxylem elements 
was not altered, whereas the number of pith cells surrounding 
the late metaxylem was significantly reduced while their size in a 
radial direction was significantly enlarged compared with cells in 
the corresponding region of wild-type primary roots (Table 1). 

Auxin control of vascular differentiation 

In a previous study, it was demonstrated that polar auxin trans- 
port was decreased by 83% in 3-d-old ruml primary roots (Woll 
et al. , 2005). To test the effect of auxin on vascular development 
in maize primary roots, wild-type and ruml seedlings were ger- 
minated either in the presence of the auxin transport inhibitor 
NPA (10 uM), the synthetic auxin aNAA (0.1 uM) or water. 
Primary roots of 60 mm length were sectioned in 10 mm incre- 
ments (Fig. 4). In the mutant ruml, disorganized vascular pat- 
terns and enlarged pith cells were observed irrespective of water, 
NPA or aNAA treatment (Fig. 4A / -R / ). Hence, treatment 
with aNAA, which can passively diffuse into the cells, did not 
recover the wild-type vascular phenotype in ruml primary roots 
(Fig. 4M / -R / ). By contrast, while aNAA treatment of wild-type 
primary roots did not alter the vascular system (Fig. 4M-R), 
NPA treatment led to a similar disorganization in vascular pat- 
terning and enlarged pith cells surrounding the late metaxylem 
(Fig. 4J-L) as observed in the mutant ruml. This suggests that 
these defects are a consequence of reduced polar auxin transport. 

The ruml mutant displays enhanced lignification of 
aberrant pith cells 

Mapman analyses of genes differentially expressed between 
wild-type and ruml primary roots revealed that three 



genes which encode key enzymes of lignin biosynthesis 
(4-COUMARATE-CoA LIGASE, 4CL; GRMZM2G16 
5844, FERULATE 5-HYDROXYLASE, F5H; GRMZM2G 
100158andCINNAMYLALCOHOLDEHYDROGENASE, 
CAD; GRMZM2G443445) were strongly induced in 2 cm 
ruml primary roots (see Supplementary Table S3 at JXB 
online). 

Secondary cell walls accumulate lignin or other secondary 
metabolites. Lignin deposition in the vasculature of wild-type 
(Fig. 5A) and ruml (Fig. 5B) primary roots was surveyed by 
phloroglucinol-HCl staining and identified by red staining 
of protoxylem and early metaxylem, and endodermis cells 
forming the Casparian strip. In contrast to wild-type primary 
roots (Fig. 5A), strong staining was detected in the larger 
and thickened pith cells surrounding the late metaxylem in 
ruml (Fig. 5B). Pith cell walls were significantly thicker in 
ruml than in wild-type primary roots (Fig. 5C), while cell 
wall thickness was not affected in ruml endodermis cells 
(Fig. 5D), which also deposit lignin. Hence, enhanced lignin 
deposition in mutant primary roots is probably not directly 
controlled by ruml but rather is a result of the defects during 
the differentiation of vascular cells. 

Discussion 

In the present RNA-Seq study, 22 833 genes were expressed 
in the differentiation zone of 2 cm primary roots of maize. 
A similar number of genes was expressed in a transcriptome 
analysis of maize primary roots by SAGE where it was extrap- 
olated from 14 850 expressed genes that ~22 000 genes are 
active in the root tips of the cultivar FR697 (Poroyko et al. , 
2005). These results suggest that at least ~55% of all high 
confidence gene models (ZmB73_RefGen_v2) are expressed 
in the differentiation zone and in root tips of young maize 
primary roots. 

RUMI controls expression of genes associated with 
auxin signal transduction 

Aux/IAA proteins such as RUMI (von Behrens et al., 2011) 
are auxin signal transduction regulators that control diverse 
aspects of plant development. In the present RNA-Seq study, 
ruml expression was significantly reduced in the differentia- 
tion zone of ruml primary roots compared with the wild type, 
confirming previous qRT-PCR experiments (von Behrens 
et al. , 201 1). Several genes known to act downstream of Auxl 



Table 1 . Quantification of xylem cell number and radial size in 80 mm WT versus rum 1 roots at a distance of 70mm from the root tip 





Number of metaxylem cells 




Pith cells around late metaxylem 






Early 


Late 


Number 


Size (um) 


WT 


13+1 


6 + 1 


110 + 7 


9.4 + 2.6 


ruml 


13 + 2 


6 + 1 


76+18 a 


20.6 + 12.5" 



a Number of WT vs ruml pith cells: p <0.05. 
b Size of WT vs ruml pith cells: p <0.01 . 
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Fig. 4. Effect of NPA and aNAA application on the formation of the vascular system of 60 mm wild-type and ruml primary roots. Series of transverse 
sections displaying the central cylinder of wild-type (A-R) and rum 1 (A'-R') primary roots grown in water (wild type: A-F; rum 1 : A'-F'), 1 0 jiM NPA (wild 
type: G-L; rumV. G'-L') and 0.1 jiM aNAA (wild type: M-R; rumV. M'-R'). Sections were taken every 10mm and show a single, representative primary 
root per genotype and growth condition. The distances are indicated with reference to the root tip. Scale bar: 50 |im. (This figure is available in colour at 
JXB online.) 



IAA genes displayed RUM 1 -dependent expression and were 
auxin-inducible. A putative role of these genes in maize root 
development is suggested by the function of their orthologues 
in root formation in other plant species. 

First, nacl (GRMZM2G081930), which encodes the 
NAC domain-containing protein 21/22 is down-regulated in 
the differentiation zone of ruml primary roots. The closest 
Arabidopsis homologue of this gene encodes the transcrip- 
tion factor NAC1 which controls lateral root initiation via 
an Aux/IAA-ARF-dependent auxin signalling module (Xie 
et al, 2000, 2002). 

Three maize plethora {pit) genes were down-regulated 
in ruml including pltl, (GRMZM2G141638), bbml 
(GRMZM2G366434), and hscfl (GRMZM2G1 39082). These 
genes encode members of the AP2/EREBP (APETALA2/ 
ETHYLENE-RESPONSIVE ELEMENT BINDING 
PROTEIN) class transcription factor family. In Arabidopsis, 
several PUT genes are down-regulated in the mutant slrliaal4 
(Aida etal, 2004) which is blocked in lateral root formation. The 
maize pltl and bbml genes are orthologues of AtBBM, which 
controls adventitious root formation (Srinivasan et al, 2007). 



Moreover, maize hscfl is the closest relative of Arabidopsis 
PLT3 and PLT7 (see Supplementary Fig. SI at JXB online), 
which are involved in the regulation of the auxin efflux trans- 
porter PIN1 (Prasad et al., 2011). Although overall pinl levels 
were not altered in the ruml mutant (Woll et al. , 2005) cell- type 
specific alterations of pinl expression might explain the reduc- 
tion in polar auxin transport by 83% in the ruml primary root. 

Furthermore, in the present study the maize like-auxl 
(laxl, GRMZM2G129413) and lax2 (GRMZM2G 149481) 
which encode auxin influx carriers, were down-regulated in 
ruml. In Arabidopsis, the AUXILAX gene family members 
AUX1 and LAX3 act concomitantly in lateral root initiation 
(Marchant et al, 2002) and emergence (Swarup et al, 2008). 

Finally, maize lbd24 (GRMZM2G075499) was down-reg- 
ulated more than 4-fold in ruml primary roots. The closest 
homologue of maize lbd24 is Arabidopsis LBD16 (Majer and 
Hochholdinger, 2011). AtLBD16 controls asymmetric divi- 
sion of lateral root founder cells in Arabidopsis (Goh et al, 
2012), which is directly regulated by the SLR/IAA14-ARF7/ 
ARF19-dependent auxin signalling module which controls 
Arabidopsis lateral root initiation (Okushima et al, 2007). 
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Fig. 5. Analyses of lignification in wild-type and rum 1 primary roots. Cross-sections of the central region of 1 0-d-old wild-type (A) and rum 1 (B) primary 
roots. Phloroglucinol-HCI-stained lignified cell walls were stained in red. Lignin was detected in protoxylem and early metaxylem elements (black arrow), 
in enlarged pith cells (black arrowhead) surrounding the late metaxylem elements (x) in ruml and in endodermis cells forming Casparian strips (blue 
arrow). (C) Pith cells surrounding late metaxylem elements in the mutant ruml display significantly thicker cell walls than their wild-type counterparts 
(n=30, ***p <0.001). (D) Cell wall thickness of endodermis cells (n=18) is not significantly different between wild-type and ruml primary roots. (A,B) scale 
bar: 1 00 jam. (This figure is available in colour at JXB online.) 



In summary, in the present study the expression of several 
genes involved in auxin signal transduction was shown to 
be RUM 1 -dependent. In other species, orthologues of sev- 
eral of these genes are known to function in root develop- 
ment. Hence, genetic dissection of these candidate genes may 
provide a better understanding of the molecular processes 
involved in maize root formation. 

RUM! controls auxin dependent xylem development 

The Arabidopsis MP/ARF5 protein, which controls vascular 
development (Przemeck et al, 1996; Hardtke and Berleth, 
1998) is orthologous to the paralogous maize auxin signal 
transduction regulators ARF8 and ARF37 (von Behrens 
et al. , 201 1 ; Wang et al. , 2012). In the present study, transcrip- 
tion of both arf8 and arf37 was almost completely repressed 
in the mutant ruml. These findings support the notion that 
transcription of arf genes is controlled by Aux/IAA proteins 
(Lau et al, 2011). To survey if arf8 and arf 37 repression 
affects xylem organization in ruml primary roots, a compara- 
tive histological analysis of ruml versus wild-type vasculature 
was undertaken. These experiments led to the observation 
that the fate of the vascular system in ruml primary roots 
depended on its developmental status. 

The vascular system of the mutant ruml formed soon 
after germination at a distance of 10 mm from the root apex 



displayed only minor defects in xylem organization (Fig. 3A', 
C, G', NT). It could be hypothesized that, at this early stage 
of primary root development, the lack of ruml expression 
is compensated by another member of the maize Aux/IAA 
gene family. This notion is supported by a recent expression 
analysis of the maize Aux/IAA gene family reporting very 
low gene expression levels of ruml but much higher levels of 
other Aux/IAA genes such as IAA14 in the meristematic zone 
of maize primary roots (Ludwig et al, 2013). Moreover, at 
this developmental stage, auxin synthesized in the root apex 
(Feldmann, 1980) could compensate for the significantly 
reduced auxin amounts transported from the shoot to the 
root (Woll et al, 2005). Subsequently formed xylem elements 
at a distance of 20^10 mm from the root apex were in severe 
disarray and their surrounding pith parenchyma cells were 
significantly enlarged (Fig. 3B', D'-F', H'-J'; N'-O'). This 
observation suggests that RUM1 controls the development 
of distinct cells-types such as xylem elements and pith cell 
parenchyma. In Arabidopsis, it was reported that a mutu- 
ally inhibitory feedback loop between auxin and cytokinin 
sets up the boundary between xylem and neighbouring pro- 
cambial parenchyma cells (Bishopp et al, 2011). Significant 
alterations in auxin and cytokinin signalling can affect the 
positions of these domains. Since the enlarged pith cells 
around the metaxylem elements are lignified, such as xylem 
cells (Fig. 5B), one could assume that these pith cells have 
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altered their cell identity. However, in Arabidopsis, protox- 
ylem cells are denned by high auxin signalling while the sur- 
rounding procambial cells are characterized by low auxin and 
high cytokinin signalling (Bishopp et al. , 201 1). For instance, 
Arabidopsis mutants with significantly impaired auxin sig- 
nalling such as axr3 do not form xylem cells at all (Leyser 
et al, 1996). Therefore, if one assumes a similar mechanism 
in maize, one would rather expect a reduction of the xylem 
elements in a mutant with compromised auxin signalling 
such as in ruml. Remarkably, the phenotype of disorgan- 
ized xylem elements and enlarged lignified surrounding cells 
is only transiently observed during a defined developmental 
period. In continuously proliferating older primary roots of 
60 or 80 mm in length, no xylem elements were formed later 
in development (60 mm: Fig. 3K'-L'; 80 mm: Fig 3P-T). 
Hence, at this developmental stage, ruml roots display the 
same low auxin signalling phenotype as observed in the 
Arabidopsis axr3 mutant (Leyser et al, 1996). 

It was demonstrated previously that the AuxlIAA mutant 
ruml displayed significantly reduced polar auxin transport 
(Woll et al, 2005). A link between reduced polar auxin trans- 
port in ruml roots and the observed defects in vasculature 
formation was demonstrated by a treatment of wild-type 
roots with the polar auxin transport inhibitor NPA. NPA- 
treated wild-type roots showed the similar disarray in vascular 
organization as observed in the mutant ruml. This supports 
the notion that polar auxin transport is important for the 
differentiation of xylem elements and the surrounding pith 
cells in maize (Jansen et al, 2012). By contrast, aNAA was 
not able to restore the wild-type phenotype in ruml primary 
roots. This is in line with the earlier observation that auxin 
treatment cannot induce lateral roots in the mutant ruml 
(Woll et al. , 2005), and supports the molecular model explain- 
ing RUM1 function (von Behrens et al, 2011). According to 
this model, the ruml mutant protein constitutively represses 
downstream auxin responsive genes even at high auxin levels 
(von Behrens et al, 2011). The same effect and phenotype 
can be obtained in wild-type plants treated with NPA which 
inhibits polar auxin transport and thereby reduces cellular 
auxin levels. 

A similar defect in vascular development has been demon- 
strated in tomato IAA15 (Deng et al, 2012). However, while 
ruml specifically controls vascular development in roots but 
not in the shoot, tomato IAA15 regulates xylem development 
in the stem (Deng et al. ,2012). 

Excessive lignification of pith cell-walls in the primary 
root of the mutant ruml 

The hydrophobic cell wall polymer lignin is deposited in 
endodermis and pith cells surrounding xylem elements to 
make them impermeable to water and to provide structural 
support (Vermerris et al, 2010). Histochemical staining 
revealed excessive lignin deposition in enlarged pith cells 
of the mutant ruml while lignin deposition in endodermis 
cells was not significantly affected. Comparative RNA-Seq 
experiments of young maize primary roots revealed strong 
up-regulation of three lignin biosynthesis genes in mutant 



ruml roots. Recent studies revealed complex transcriptional 
networks that control lignin biosynthesis in addition to a 
plethora of developmental and environmental cues includ- 
ing various types of stress and auxin (Zhao and Dixon, 
2011). The maize mutant Irtl (lateralrootlessl) displays 
a similar phenotype as ruml (Hochholdinger and Feix, 
1998). Both mutants fail to initiate lateral roots from the 
primary root. The Irtl gene has not yet been cloned but it 
has been demonstrated that it is not allelic with ruml (Woll 
et al, 2005). A comparative proteome analysis of 9-d-old 
wild-type and Irtl primary roots revealed four proteins 
involved in lignin biosynthesis (including Caffeoyl-CoA-3- 
O-methyltransferase, CCoA-OMT) that were up-regulated 
in Irtl primary roots compared with wild-type primary 
roots (Hochholdinger et al, 2004b). Similarly, a compara- 
tive transcriptome analysis of wild-type versus ruml pericy- 
cle cells revealed preferential expression of CCoA-OMT in 
mutant pericycle cells (Woll et al, 2005). 

Moreover, it has been demonstrated that inhibition of 
auxin transport in Arabidopsis resulted in dramatically 
increased numbers of cells with thickened secondary cell 
walls in stems which suggests that cell wall thickness is also 
regulated by auxin (Mattsson et al, 1999). Therefore, drasti- 
cally decreased polar auxin transport in ruml might lead to 
the mis-regulation of lignin biosynthesis genes resulting in 
thickened secondary cell walls in enlarged pith cells. 

In summary, the data presented here indicate that the tran- 
scriptional repressor RUM1, in addition to its previously 
demonstrated function in regulating lateral and seminal root 
initiation, also controls vascular and pith cell differentiation 
and lignin deposition in these cells. 

Supplementary data 

Supplementary data can be found at JXB online. 

Supplementary Fig. SI. Phylogenetic tree of nine maize 
and eight Arabidopsis pit genes generated by ClustalV. 

Supplementary Table SI. Summary of the RNA-Seq data 
of ruml and the wild type, and alignments to the B73 refer- 
ence genome (ZmB73_RefGen_v2). 

Supplementary Table S2. List of the 22 833 expressed genes 
and their characteristics. 

Supplementary Table S3. Differentially expressed genes 
(Fc >2; FDR <1%) were functionally annotated using the 
MapMan software. Genes that were assigned to more than 
one functional class or subgroup are labelled in red. Fold 
changes (ruml/WT) are given as logarithmic (log 2 ) values. 

Supplementary Table S4. Determination of overrepre- 
sented and underrepresented functional classes among differ- 
entially expressed genes. 

Supplementary Table S5. Sequences of oligonucleotide 
primers used for qRT-PCR analyses. 
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